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ABSTRACT 



I . Context. Magnetars are an extreme, highly magnetized class of isolated neutron stars whose large X-ray luminosity is believed to be driven by 
. ■ their high magnetic field. 

t-H ' Aims. Study for the first time the possible very high energy y-ray emission above 100 GeV from magnetars, observing the sources 4U 0142+61 
y—i ', and IE 2259+586. 

O^l ■ Methods. We observed the two sources with atmospheric Cherenkov telescopes in the very high energy range (E> 100 GeV). 4U 0142+61 was 

t— I ■ observed with the MAGIC I telescope in 2008 for -25 h and IE 2259+586 was observed with the MAGIC stereoscopic system in 2010 for ~14h. 

!LJ ' The data were analyzed with the standard MAGIC analysis software. 

, , Results. Neither magnetar was detected. Upper limits to the differential and integral flux above 200 GeV were computed using the Rolke algorithm. 

. We obtain integral upper limits to the flux of 1.52xlO~' 2 cirT 2 s and 2.7xlO~ 12 cirT 2 s with a confidence level of 95% for 4U 0142+61 and 

% , • IE 2259+586, respectively. The resulting differential upper limits are presented together with X-ray data and upper limits in the GeV energy range. 
2 ■ 

. . . ■ Key words. Radiation mechanisms: non-thermal, Stars: magnetars, Gamma rays: stars, Stars: individual: 4U 0142+61, IE 2259+586 



1. Introduction 

Magnetars are a peculiar class of neutron stars. Most of the about 
20 known magnetars are characterized by strong dipolar mag- 
netic fields (~ 10 14 - 10 15 Gauss) that are ~ 10 - 1000 times 
higher than the average value in radio pulsars, near or even above 
the quantum electrodynamic field strength, Bqed = m^c 3 /efi ~ 
4.4 x 10 13 G ((Harding & Lail 120061) . although with two excep- 
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tions dRea et ai]|2010l |2012|) . They have bright X-ray luminosi- 
ties (L x ~ 10 32 - 10 36 ergs~ 1 ) from 0. 1-300 keV, longer rota- 
tion periods than most ordinary radio pulsars (~ 2-12 s), and 
very high period derivatives (~ 10~ 13 -1 0~ n s/s). For mor e de- 
tails see recent review s on magnetars by Mereghetti (2008f) and 
lRea& Espositol (120111) . 

The most successful model for explaining the X-ray 
emission from these objects in vokes the decay and insta- 
bility of their magnetic fi e lds (Ehincan & Thompson! | l992; 
Thompso n & Duncan! 1 19931 1 19951) . The dichotomy between 
magnetars and ordinary pulsars may indicate different progen- 
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itors (iThompson & Duncad fl996). This scenario for birthing a 
magnetar postulates a very rapidly spinning proto-neutron star at 
birth, which would then have a high rotational energy. This ex- 
cessive energy was searched in the supernova remnants (SNRs) 
surrounding magnetars. However, those SNRs sh ow no excess in 
X-ray s relative to those around normal pulsars dVink & Kuipeil 
120061) . 

With this paper we aimed at testing whether the additional 
energy at birth could have gone in TeV emission. 

To date, no magnetar has been detected at energies above 
1 MeV. Although recently, the H.E.S.S. collaboration presented 
their discovery of extended TeV y-ray emission towards the 
magnetar SGR 1806-2CQ, it i s doubtful that the emission is 
driven by the magnetar itself dRowell et alj|201 lb . The Fermi- 
LAT Collaboration presented upper limits for 13 magnetars after 
17 months of sky s u rvey observations betw e en 0. 1 and lOOGeV 
dAbdo et al.ll2010l) . ISasmaz Mus & Gogiisl d2010t) studied spe- 
cially the Fermi data of 4U 01 42 +61. Neither steady nor pulsed 
emission was found. In this work we present a search for the 
emission at very high energies (200 GeV-50 TeV) from the two 
magnetars 4U0142+61 and 1E2259+586 with the MAGIC tele- 
scopes. These sources have been also observed by the VERITAS 
Collaboration and corresponding u pper limits above an e nergy 
of 400 GeV have been presented in lGuenette et all (120091) . The 
present MAGIC observations of these two magnetars extend the 
spectrum to lower energies, 200GeV. 

2. The observed magnetars 

The source 4U 0142+61 is located at Q^ooO; ^2000 — 
01 h 46 m 22H07,+61°45'03'/19 at a distance of ~ 3.6 kpc. 
With an X-ray luminosity of Lx ~ 1 x 10 35 erg s -1 it is one of the 
most X-ray luminous magnetars known (Mc Gill Pulsar Group! 
120121) . This makes it a good target to search for persistent 
very high energy emission. Long term spin period variations 
P ~ 8.7 s) were d iscovered during observations with EXOSAT 
Israe l et all 1 1994b . leading to the measurement of the period 
derivative P ~ 2 x 10~ 12 ss~', and consequen tly of the very 
strong magnetic field B ~ 1.3 x 10 14 G dMcGill Pul sar Group 
12012b . The bright l-10keV emission coming from 4U 0142+61 
has been observed by many X-ray satellites dWhite et al.l l 1987; 
llsrael et alfl999t iPatel et alj|2003l: iRea et al.ll2007albl) revealing 
an X-ray spectrum typical of an Anomalous X-ray Pulsar 
(AXP), best described by an absorbed blackbody plus a power 
law (N H ~ 10 22 ctrT 2 , kT ~ 0.4keV and T ~ 3.62). A very strong 
hard X-ray emission has been reported by INTEGRAL up to 
250 keV, with a spectrum wel l modeled wi t h a st eep power-law 
with a photon index of ~1 dKuiper et al.l l2006h . At the time 
of data taking with the MAGIC telescope, there were only 
COMPTEL upper limits in the MeV range suggesting a spectral 
break in the hard X-ray emission of this object. The upper 
limits, however, do not put strong constraints on the HE or VHE 
gamma-ray emission of the object, especially given the high 
systematic uncertainty of the background subtraction in the data 
COMPTEL analysis dSchonfelderil2004l) . Recently the upper 
limits derive d by the Fermi-LAT Collab oration dAbdo et alj 
120101) and bv lSasmaz Mus & Gog~i is (2010) point to a cutoff in 
the MeV band. 

The AXP IE 2259+586 is located at ff 2 ooo, £2000 = 
23 h 01 m 08 s .296, +58°52'44'/45 embedded in the SNR CTB109. 



1 Originally, magnetars were divided into two categories: Soft 
Gamma Repeaters (SGRs ) and Anomalous X-ray Pulsars (AXPs) 
dWoods & ThompsorJl2006n . 



The source has a magnetic field of B ~ 0.59 x 10 G and a 
distance of ~4 kpc, making it a good c andidate for MAGIC ob- 
servations dMcGill Pulsar Gro up 2012). RXTE measured the pe- 
riod (P ~ 7 s) and the period derivative (P ~ 0.5 x 10~ I2 ss~') 
dGavriil & Kaspi 2002). The X-ra y spectrum is var i able depend- 
ing on the source emission state (Kaspi et al. 2003t lWoods et al.l 
2004). After undergoing an outburst in 2002, the source re- 
turned into its possible quiescence state and the correspond- 
ing spectrum is best fitted by a blackbody plus a power la w 
(N H ~ 10 22 cnr 2 , kT ~ 0.4keV and F ~ 3.75) (IZhu et al]l2008h . 

3. The MAGIC telescopes, analysis and data 

The MAGIC Collaboration operates two 17 m diameter imag- 
ing atmospheric Cherenkov telescopes on the Canary Island of 
La Palma. The data sets presented here were taken in 2008, 
i.e. before the second MAGIC telescope was operational (mono 
data), and in 2010 when both telescopes were already taking 
stereoscopic data. Details about the perfo rmance of M AGIC 
in mono and stereo mode can be found in lAlbert et al.l (120081) 
and lAleksic et al.1 d2012l) . All data presented in this work were 
taken in the so-called wobble mod e and were analyzed using the 
MARS analysis framework (Mor aleio et alj|2009l lAleksic et al.l 
2012). The analyses presented here have an analysis threshold 
of 200 Ge V. The upper Umit s were calculated using the Rolke 
algorithm dRolke et al J 12005) with a confidence level (C.L.) of 
95% assuming a Gaussian background and 30% of systematic 
uncertainty in the flux level. Since IE 2259+586 is embedded 
in a SNR and may contain more than one emission region (see 
below) relevant parameters for the observations are the MAGIC 
field of vi ew of 3.5° and the a ngular resolution of ~0.07° above 
300 GeV (lAleksic et al.ll20H . 

4U 0142+61 was observed for 25.41 hours. After quality cuts 
16.58 hours of effective observation time remain. These mono 
data were taken between August and December 2008 covering a 
zenith angle range between 33° and 40.6°. 

Data for IE 2259+586 were taken in stereo mode wobbling 
around the sky position 0.12° away from the magnetar to have 
the shell of the supernova remnant and the magnetar in the same 
field of view. Given the angular resolution of the MAGIC tele- 
scopes, these two possible TeV sources would be spatially sep- 
arable with MAGIC. The region was observed between August 
and November 2010 for 14.33 hours within a zenith angle range 
of 29°-43°. After quality cuts this amounted to 8.22 hours of ef- 
fective observation time. 

4. Results 

Neither source was detected by MAGIC. We computed the inte- 
gral flux upper limits above 200 GeV with 95% C.L. assuming a 
differential energy spectral shape of a power law with an index 
of 2.6, similar to that of the Crab Nebula spectrum. The results 
are given in Table [TJ A 25% change in the photon index yields 
a variation of about 7%. In Fig. [T]we show the corresponding 
test statistic (TSfl map for IE 2259+586. No excess was found 
at either the magnetar position nor at any location within the 
surrounding SNR. The TS map for 4U 0142+61 is not shown 
here, but shows the same flat behaviour. The upper limit for the 
extended SNR will be discussed elsewhere. The white contours 
represent the X-ray emission of the surrounding SNR detected 

Our test statistic is iLi & Mai [1983) eq. 17, applied on a smoothed 
and modeled background estimation. Its null hypothesis distribution 
mostly resembles a Gaussian function, but in general can have a some- 
what different shape or width. 
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Table 1. Magnetar parameters taken from McG ill Pulsar Groupl (120121) . along w ith the MAGIC resu lts presented here. Crab Units 
(C.U.) are defined as a fraction of the Crab Nebula flux as measured by MAGIC (Aleksic et al. 2 0T2h . 



Source 


Distance 


Ssurf 


Lx 


log(Aot) 


Eff. obs. time 


Significance 


Upper limit (95% C.L., E>200GeV) 




[kpc] 


[10 14 G] 


[10 35 erg s- 1 ] 


[erg s _1 ] 


[hrs] 


<x 


[cirr 2 s -1 ] 


4U 0142+61 


3.6+0.4 


1.3 


1.1 


32.10 


16.58 


-2.1 


1.52X10 -12 (0.70% C.U.) 


IE 2259+586 


4.0+0.8 


0.59 


0.34 


31.70 


8.22 


-0.5 


2.70xl0~ 12 (1.24% C.U.) 



with the XMM-Newton satellite (0.1-15keV). We also searched 
for pulsations for both magnetars. For the pulsed analysis of IE 
2259+586 we used a timing soluti on valid at the epoc h of the 
MAGIC observations, as derived bv llcdem et alj(l2012l) . We did 
not detect any significant pulsation at VHE energies. In the case 
of 4U 0142+61, we searched fo r pulsation using the ephemeris 
of ISasmaz Mus & Gogiisl d2010h . We did not find any pulsation 
at VHE energies for this source either. 

Since neither source experienced an outburst in X-rays 
during our observing intervals, we can compare our up- 
per limits with data taken with different instruments during 
different quiescent epochs. In Fig. [2^ (HJ3) we present the 
0.1keV-3TeV multi-band spectral energy distribution (SED) 
of 4U0142+61 (1E2259+586), respectively. For both sources 
the corresponding differential and integral upper limits derived 
in this work are shown (red lines in Fig. [2). In the case of 
4U 0142+61, the 0. 1-200 k eV data are from XMM-Newton- 
PN a n d INTEGRAL-ISGR1 dRea et all l2007at Iden Hartog etail 
2008; iGonzalez etail 120101) plotted toge ther with the 2cr 
COM PTEL upper limits dden Hartog et al I l2006t iKuiper et alj 
120061) . For IE 2259+586 we s how data points from XMM- 
Newton-P N dWoods et al.ll2004l) together with COMPTEL up- 
per limits ( Kuiper et alj I2006I) . The upper limits provided by 
the Fermi-LAT C ollaboration were calculated for three differ- 
ent energy ranges ( Abd o et ai1l20 1 Oh . For the overall energy bin 
from 0.1-10GeV a photon index of 2.5 was assumed. A cut- 
off is mimicked by splitting this energy bin into two parts with 
photon indices of 1.5 and 3.5, respectively. The assumed slopes 
are indicated in Fig. [2] The results derived by the VERITAS 
Collaboration on the two sources are also shown for compari- 
son (grey dashed lines). They correspond to 99% C.L. integral 
flux upper limits of 8.68 x 10~ 13 cnT 2 s _1 for 4U 0142+61 and 
2.49 x 10 _12 crrr 2 s _1 for IE 2259+586 by assu ming a power- 
law w ith a photon index of 2.5 above 400 GeV dGuenette et alJ 
2009). The upper limits for both sources are compatible with a 
break in the power law at ~1 MeV. However, the SED lacks any 
measurements above hard X-rays, what gives complete freedom 
under the corresponding instrumental sensitivity. 



ICheng & Zhangl (l200lh presented a model for the very high 
energy radiation from magnetars. They predicted emission of 
y-rays in the GeV band coming from the outer gap for the 
tw o sources we stu died. This model has been recently revised 
by iTong et ail d201 lb . who updated the observational param- 
eters to calculate the y-ray radi ation properties of all A XPs 
and SGRs using the models by IZhang & Chengl (Il997h and 
ICheng & Zhang dJOQl). The scenario bv lTong et al.l (|201 1) pre- 
dicts that4U0142+61 should have been detected by ferm/-LAT, 
although they explain the lack of a detection by Fermi-LAT 
(Abdo et aI] |2010l ; ISasmaz Mus & Gogtis 2010) by invoking ac- 
cretion. For IE 2259+586 the model does not predict GeV emis- 
sion. We note that although none of the current models predict 
TeV range emission for either magnetar, the existence of dif- 
fuse emission around IE 2259+586 could lead to the appearance 
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Fig. 1. TS map of 1E2259+586. The green cross represents the 
magnetar position. The white contours show the X-ray emis- 
sion of the surrounding SNR CTB 109 detected with the XMM- 
Newton satellite. The color scale represents the TS value. 



of an extra component in the SED besides any magnetospheric 
emission. 

Using the MAGIC telescopes we studied for the first time the 
possibility of magnetars to be a new TeV source class on the ex- 
amples of 4U 0142+61 and 1E2259+586. This exploratory work 
led to a non-detection of the VHE gamma-ray emission from 
either of them. This result indicates that magnetars are proba- 
bly not VHE emitters during their quiescent state, as expected 
from the various theoretical models. However, the possibility of 
magnetars being VHE emitters during flaring episodes cannot 
be ruled out because of the lack of VHE observations during 
these high-activity periods. Consequently, our future searches 
for VHE emission of magnetars will be performed during out- 
bursts 0. 

Acknowledgements. We would like to thank the Institute) de Astrofi'sica de 
Canarias for the excellent working conditions at the Observatorio del Roque 
de los Muchachos in La Palma. The support of the German BMBF and MPG, 
the Italian INFN, the Swiss National Fund SNF, and the Spanish MICINN is 
gratefully acknowledged. This work was also supported by the CPAN CSD2007- 
00042 and MultiDark CSD2009-00064 projects of the Spanish Consolider- 
Ingenio 2010 programme, by grant DO02-353 of the Bulgarian NSF, by grant 
127740 of the Academy of Finland, by the DFG Cluster of Excellence "Origin 



3 In order to provide fast reactions to such events in the future, 
MAGIC has installed an alert system, which receives alerts provided 
by several satellites and points the telescopes to the flaring source auto- 
matically, as it is also done for observations of Gamma Ray Bursts. 
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